Theory of Magnetocrystalline Anisotropy Energy for Wires and Corrals of Fe 

adatoms: A Non-Perturbative Theory 
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The magnetocrystalline anisotropy energy E an is for free-standing chains (quantum wires) and 
rings (quantum corrals) of Fe-adatoms N =(2. ..48) is determined using an electronic tight-binding 
theory. Treating spin-orbit coupling non-perturbatively, we analyze the relationship between the 
electronic structure of the Fe d-electrons and E an is(rid), for both the chain and ring conformations. 
We find that Eanis (N) is larger for wires than for rings or infinite monolayers. Generally E a nis(nd) 
decreases in chains upon increasing N, while for rings E an is(nd) is essentially independent of N. 
For increasing N, E a ni s {nd) in corrals approaches the results for freestanding monolayers. Small 
rings exhibit clear odd-even oscillations of E an i S (N). Within our theoretical framework we are able 
to explain the experimentally observed oscillations of E an i s (rid) during film growth with a period 
of one monolayer. Finally, a generalization of Hund's third rule on spin-orbit coupling to itinerant 
ferromagnets is proposed. 
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I. INTRODUCTION 



Over the past decades information technology has been 
determined by semiconductor applications. However, re- 
cently the technology of such devices is limited due to 
the optical wavelenght. Thus, quantum structures on the 
nanoscale can offer new perspectives for electronic appli- 
cations. Manipulations at this scale require to develop 
nanotools: Low-temperature scanning tunneling micro- 
scopes allow not only for the probing of metallic surfaces 
with atomic resolution, but also for single-atom manip- 
ulation on surfaces. It has recently been demonstrated 
that lateral structures of adatoms can be generated by 
the scanning tunneling, microscope (SjTM) tip (Crommie 
et aim, Meyer et alB, Rieder et alE). Modern STM- 
techniques allow now for the preparation of quantum cor- 
rals (QC's). Examples of QC's which have been prepared 
are rings of 48 Fe adatoms U or stadia El. The latter shape 
is motivated by the search for quantum-chaotic behav- 
ior, which so far has not been successful due to strong 
scattering of the electrons by the corral atoms. However 
the magnetic moment and magnetocrystalline anisotropy 
for these nanoscopic structures have not yet been deter- 
mined in experiment. The QC of adsorbate atoms sets 
the boundary condition for the wave function of the sub- 
strate surface state. The quantum mechanical electron 
distribution of this state can therefore be measured by 
STM and is found to be in remarkable agreement with 
standard predictions of quantum mechanics. 

In this paper, however, we are not interested in the 
surface state but rather in the magnetic anisotropy prop- 
erties of the quantum corral itself. Since it is known 
that the magnetocrystalline anisotropy energy E an i S is 
enhanced in thin magnetic films compared to threedimen- 
sional bulk media, one might speculate about a further 
increase of E an i S upon reduction of the dimensionality 



to ID. Thus, we calculate E an i S for QC's and quantum 
wires (QW's). These chains are frequently realized by 
the enhanced adsorption probability of adatoms on step 
edges. 

We analyze for chains and rings for different numbers 
N of Fe-adatoms how the magnetocrystalline anisotropy 
energy E an i S depends on the electronic structure. To our 
knowledge, the -first attempt of this kind has been made 
by Wang et alu for a pair of two Fe-atoms to explain 
the dependence of the electronic states. While Wang et 
al. treated spin-orbit coupling (SOC) in second-order 
perturbation theory, we include SOC completely non- 
perturbatively without resorting to degenerate or non- 
degenerate perturbation theory of any order. It is shown 
that for the diatomic pair there are important contri- 
butions to the magnetocrystalline anisotropy energy as 
a function of the number nd of d-electrons per atom 
E an is{nd) due to the different lifting of band degenera- 
cies for different magnetization directions . We find that 
Eanis oc Ag Q is valid for the diatomic pair. However, 
E a nis(n>d) can change its sign with respect to A so due 
to level crossings in the energy level scheme of the di- 
atomic pair. In the same way we calculate E an i S (nd) for 
chains and rings with different numbers N of Fe-atoms. 
In small rings there are oscillations of E an i S (N) between 
rings with an even or odd number of atoms. For our cal- 
culation we use two models: n 

(i) an extension of the diatomic pair modeld to N atoms, 

(ii) a conventional tight-binding description of the wires 
and corrals including SOC. 

Our paper is organized as follows: In Sec. II, the 
diatomic-pair model (II. A) and the tight-binding model 
(II. B.) are presented together with the determination of 
SOC (II. C) and E anis (II. D.). Results for the diatomic- 
pair model are shown in Sec. III. A, while results for the 
tight-binding model are discussed in Sec. III.B. 
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II. THEORY 



A. Diatomic-pair model 



We follow the ideas of Wang et al. and extend their 
diatomic pair model to the case.-pf N atoms. We start 
from the following HamiltonianEfQ: 



^ilj^ C t 1 a C jl3a' (•*") 



The summation is taken over N sites i, while the sec- 
ond sum is over pairs of nearest neighbors. As usual 
ci ff (cj 7CT ) creates (destroys) an electron with spin a at 
site i and orbital 7. Furthermore £i 7 is the on-site energy 
of an electron at site i and orbital 7. For simplicity we 
set 6i 7 =0, since we consider only orbitals of the 3d type. 
The basis functions k = 1, 5 (k = 6, 10) are the 
atomic 3d orbitals conventionally denoted by xy, yz, zx, 
x 2 — y 2 and 3z 2 — r 2 , where z denotes the interatomic 
axis, together with the spin eigenstates ||) and ||) with 
respect to the spin quantization axis zm- The complete 
ferromagnetic tight-binding Hamiltonian is 



Htb — H + H ex , 
with the diagonal on-site matrix H, 

rW - 



ex 



ex ex 



J, 
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ex and 



(2) 



(3) 



which is actually independent of i. 1 and —1 denote 
5x5 diagonal matrices, and O is a 5 x 5 zero matrix. 
The matrix of the Hamiltonian H t b is then given by: 



H, 



n, 
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ex 
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with the intcrsite Hamiltonian matrix 
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(5) 



diagonal elements of the off-diagonal Hamiltonian flfc'Z. ■ 
Consequently the nearest neighbor hopping U^jp of Eq. 
(1) reduces to ti 7J - 7 , which is equivalent to the three in- 
dependent nonvanishing interactions Vdda, Vddir, VddS in 
Eq. (6). In rings, with periodic boundary conditions, 
the interaction between the first and last atom of the 
corresponding chain must be additionally included. As 
in ref. [5] we use V dd a = -0.25 eV, Vm* = 0.18 eV, 
VddS = —0.04 eV, J ex = 3.0 meV, and lattice parameter 
a=5.98 a.u. which matches the fictitious W(001) sub- 
strate. 



B. Tight-binding model 



The main disadvantage of the diatomic pair model is 
its unrealistic description of tilted <i-bonds. To perform 
more realistic calculations of E an i S (rid) in rings, we have 
to consider how the bond angles are changed between 
nearest neighbors of Fe-atoms. This aspect is taken into 
account in the tight-binding model, within the two-center 
approximation u. The tight-binding Hamiltonian H t b has 
the same structure for this model as the matrix in Eq. (4) . 
Only the elements of flfci will change their values and 
positions, because the bond angles have changed. The 
Hamiltonian is parametrized in terms of the two-center 
integrals Ed,d(l, m, n) of Slater and KosteiUJ, given in Ta- 
ble I, with d as one of the five 3d-functions xy, yz, zx, 
x 2 — y 2 and z 2 , the latter denoting 3z 2 — r 2 for simplic- 
ity. The direction of the vector Rj — Ri , pointing along 
the bond from one atom to each of its nearest neighbors, 
is given by the direction cosines I, m, and n. Thus in 
the tight-binding model both blocks 



= H 



\ 1 



of 



(4) H zntl change to: 
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VddS 





V 











Vdda I 



(6) 

Due to the symmetry of the atomic 3d-orbitals, the 
dominant overlaps result from the d-orbitals of the same 
kind, thus yielding nonvanishing contributions only in the 



Since we take the (x, y)-plane as the ring or chain plane, 
all direction cosines n must be zero, thus further sim- 
plify the nonvanishing two-center integrals of Table I. 
The other two direction cosines I and m are determined 
by the bond angle between the Fe-atoms. 

If we calculate the magnetic anisotropy of single or 
double chains we assume them to lie also in the (x, y)- 
plane and being parallel to the x-axis. We couple 
them either ferromagnetically or antiferromagnetically, 
by changing the sign of J ex . 
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C. Spin-orbit coupling 



D. Anisotropy energy 



Magnetocrystalline anisotropy is caused by SOC be- 
tween the d-states. SOC is introduced in the usual spher- 
ical on-site form as H so = A so l ■ s, with the orbital and 
spin moment vectors 1 and s respective and the atomic 
spin orbit coupling constant A so . Expressing the compo- 
nents of the orbital momentum operator 1 in the rotated 
frame (xm,Um, zm), where zm is the spin quantization 
axis, which is parallel to the direction of magnetization 
.we may rewrite H so = ^2 Hso in the following 
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(10) 



w, u 



', and w' are given by 
sin #m sin 0m, w = cos 9m, u' = 



t>M,v = cos <pM+i cos 6*m sin (p M , 



The variables u, v 
u = sin #jvf cos <pM , w : 
— sin 4>M+i cos 6* m cos 
and w' = — isin^M- 

The value for the SOC constant A so is taken from the 
Fe-atonM X so = 50 meV. 

Now we have set up the Hamiltonian matrix in terms 
of the magnetization direction (9,(f>). Since SOC oper- 
ates to a good approximation only on-site, the matrix 
elements of H\ Q occur exclusively in the diagonal blocks 
of H t b in Eq. (4); consequently, the matrix elements re- 
ally occur only in the diagonal and off-diagonal blocks of 
Hex ■ For chains, Hso is independent of site i, while for 
rings the situation is more complicated. If the magnetic 
moments are all parallel, Hso is independent of the site 
index i. However, if the magnetic moments have a radial 
configuration, Hso is dependent on i. 
Now the total Hamiltonian H to t is 

Htot = H + H ex + H so . 

In our treatment SOC is included non-perturbatively 0. 
So we get also contributions to E an is(rid) from higher 
order perturbation theory and may obtain information 
about the scaling behavior of E an i S (n d ) with the SOC 
constant A so . 



In both models we define the magnetic anisotropy en- 
ergy as the difference of energies for perpendicular mag- 
netization E±_ and parallel magnetization E\\ : 



E a nis(nd) '■— E± — £j| 



(11) 



Note that for the diatomic pair-model (as defined in Sec. 
2A) E± (E\\) refers to magnetization perpendicular (par- 
allel) to the interatomic bonds. For the tight-binding 
model (as defined in Sec. 2B) E± (En), however, refers 
to magnetization perpendicular (parallel) to the plane of 
the QC, which is the (x,y)-plane. The QW's in this case 
are assumed to lie in x direction. Thus E± (En ) refers to 
M || z (M || x). 

We define the magnetic anisotropy energy per atom for 
the diatomic-pair model as: 



E„ 



,(n d ) := E tot (0 = tt/2, ; n d ) - E tot (8 = 0; n d ). (12) 



Due to the chosen coordinate system, there is no in- 
dependence of E an i s (rid) in contrast to the tight-binding 
model. For the tight-binding model we define E an i S (rid) 
as: 

E a nis(nd) ■= E tot (0 = 0;n d ) - E tot (0 = 7r/2,(f>;n d ), 

(13) 

The in-plane angle (f> is choosen such that the result- 
ing 1-EqmsI is the largest possible. E to t(9, (f>, n d) is the 
ground-state energy per atom of the chain or ring with 
a total of rid 3d-electrons per atom, N Fe-atoms and the 
magnetization direction denoted by (9, <p) u. The total 
energy per atom E to t is given by 

E to t(9, <f>; n d )= ^ J2 E ™(°> CM*' $ ~ E ^ 9 ' & nd ^ 

m 

Here is /o(A£?) the Fermi-function at zero temperature 
and Ep(9,(/)]nd) is the Fermi-energy. For a given band- 
filling rid, Ep(9, (j>;rid) is determined self-consistently by 



rid 



— f (E n 



(6,<p)-E F (6,c(>; n d )), 



The m-th eigenvalue of the total Hamiltonian Htot with 
magnetization along (9, <fi) is given by E m . 

Note that the magnetic-moment direction enters only 
by SOC and can be adjusted for each atomic site indi- 
vidually. 

Concerning the anticipated magnetic order of the QW 
and QC in our electronic theory we remark the following: 
(i) The interatomic distance in our model is chosen in ac- 
cordance with the lateral lattice constant of a thin bcc- 
Fe film on a W(001) substrate, which is magnetically 
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ordered. However, the STM-experiments produce QC 
with the interatomic distance of 8-10 rA, which proba- 
bly would not order, except for possible indirect-exchange 
effects related to the substrate. Our interest, however, 
is mainly in small QC's, which might be produced in 
the near future, with closer nearest-neighbor distances 
and magnetic order. Furthermore it is important to note 
that the assumption of magnetic long-range order in our 
electronic theory is not in contrast to the well known 
general-theorems of the (ii) Ising-liij and (iii) Heisenberg- 
modelstj. 



III. RESULTS AND DISCUSSION 
A. Results from the diatomic- pair model 

In Fig. [l], the magnetocrystalline anisotropy energy 
Eanisirid) per atom is shown as a function of 3d- 
bandfilling rid for the Fe diatomic pair (dashed curve). 
First we note that, as expected, E an i S is enhanced com- 
pared to Fe monolayers (E an i S = 1 ... 16 meV for the 
dimer vs. ss 300 /xeV for an Fe monolayer), because of 
the symmetry reduction from 2D to ID. This is in accor- 
dence with ipcent studies of magnetism in nickel clusters. 
Apsel et aZ.til observe enhanced magnetic moments for 
small nickel clusters and find magnetization minima for 
clusters with closed geometrical shells and maxima for 
relatively open clusters. 

In Fig. [l], also the dependence of E anis (rid) on the scal- 
ing of all d-electron hopping parameters with a overall 
parameter t for the diatomic pair is presented (dotted 
curve). The absolute value of the magnetic anisotropy 
energy \E an i s (rid)\ increases for decreasing t (decreasing 
bandwidth) for the easy axis perpendicular to the bond 
{Eanis {nd) < 0), and decreases for the easy axis paral- 
lel to the bond {E an i S {rid) > 0). Note that, contrary 
to the situation in thin films, no convergence problems 
occur in E an i S due to our finite system. Our calcula- 
tion is performed in real space, and we sum over dis- 
crete energy levels. Thus the change of the occupa- 
tion number is And = 0.5, and we get 20 points for 
E a nis{ n d) of the diatomic pair. When the band is com- 
pletely empty or filled, no contribution to E an i S results: 
E a nis{nd = 0; 10) = 0. While this is trivially the case 
for an empty band, this relation holds with 14 significant 
digits for the completely filled case, which yields an in- 
dependent check of our numerical accuracy. We find for 
Fe, with its atomic band filling rid=6, an in plane easy 
axis for E an i S {Fe)—6.1 meV. Since the diatomic pair is 
a small system, it is not necessary to assume implicitly 
hybridization with s-electrons. 

We find for the diatomic pair a mirror symmetry of 
E an is{ n d) with respect to half filling {rid — 5), where the 
spin-up subband is completely full. However, the cor- 
responding points of the curve are not totally symmet- 
ric, and a close inspection shows that there are devia- 



tions of this mirror symmetry of the order 10 -5 eV. This 
small symmetry violation results from the combined ac- 
tion of SOC and hopping interaction. A detailed analysis 
shows that this asymmetry increases monotonously with 
increasing SOC. Especially for rid— 2-5 and 7.5 this asym- 
metry is clearly visible in Fig. |l|. To illustrate the origin 
of this asymmetry we split the spin-orbit coupling matrix 
H so into two parts H^ r and ,ff*J 1 ** por , for the coupling 
between parallel and antiparallel spin states, and recal- 
culate E an i S {rid) with either of the two matrices instead 



HJJ 
ffii 



H n 

w - L± SO 





rrpar I TTCintipa 

■ J ■*■ « n I - ^ a n 



So we obtain the curve E a J^^ ar {rid) (solid curve) of 
Fig. ^, which exhibits a mirror symmetry with respect 
to half filling {rid = 5). Due to the large exchange split- 
ting J ex , which completely separates the spin subbands, 
Eal^l Par { n d) is very small and thus E anis sa E^ r is {n d ). 
The E a l "ll par {rid) contribution prefers an easy magnetic 
direction perpendicular (parallel) to the molecule axis 
if the spin-down subband is less (more) than half-filled, 
and vice versa for the majority subband in agreement 
with the result of Wang et al. [5] . Also for larger chains 
Eanil par \ n d) changes it sign at rid— 2.5 and 7. 5_ which 
seems to be a general trend also in Fe-monolayersEEl. This 
is a generalization of Hund's third rule on SOC which ap- 
plies to itinerant magnetic systems (see Appendix). 
While for Fe parameters the relation E^- S {rid) + 

Eani7 a7 { n d) ~ E ams {n d ) holds to a good approxima- 
tion, small deviations are responsible for the asymmetry 
mentioned above. This asymmetry only occurs if we take 
into account at least the matrix elements of H so for the 
7r- and 5-bonds and is absent for all other combinations. 

The electronic origin of E anis {rid) of the diatomic pah- 
can be explained from its energy level scheme in Fig. 
for the two magnetization directions E m {9=0) (paral- 
lel to molecule axis) and E m (8—pi/2) (perpendicular 
to molecule axis). Only the spin- up states are shown 
in Fig. ^ for simplicity; due to the exchange splitting 
J ex =3.0 eV the spin-down states (not shown in the fig- 
ure) are located symmetrically to these bands with re- 
spect to the zero of the energy scale. The individual 
levels of the spin-up subband have been labeled by the 
numbers 1. . .10 in the absence of SOC. Without SOC, 
the lowest bonding (level 1) and the highest antibond- 
ing states (level 10) result from the splittings caused by 
Vdda- Due to the inequalities \V d da\ > | Vddvr | > \Vdds\ 
(Sec. II. B) the lowest bonding- and antibonding states 
are split by 2|Vdd<r|, while in each case two doubly degen- 
erate 7r-(5-)bonding and antibonding states result from 
Vdd-K and VddS (see Eq. 6). Thus the levels 2 and 
3 (4 and 5) of Fig. || correspond to the 7r-((5-)bonding 
states, while the levels 8 and 9 (6 and 7) are equiva- 
lent to the corresponding antibonding states. These re- 
maining degeneracies are lifted very strongly by SOC for 
E m {0=0) and therefore level crossings may occur. For 
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E m (0=pi/2) the degeneracies are only weakly lifted. The 
electronic origin of -Bonis i n d) for the diatomic pair results 
from (i) different SOC-induced shifting of occupied, 
nondegenerate levels for the two magnetization direc- 
tions, and (ii) different lifting of degenerate levels. 
Note that the lifting of degeneracies can favor a mag- 
netization parallel as well as perpendicular to the bond 
axis. 

Next we investigate the dependence of E anis on X so . 
We find E anis (nd) oc A^, for all of the nd values of the 
diatomic pair in agreement with Wang et alB, irrespec- 
tive of the value of the spin-orbit coupling constant. To 
demonstrate this, we show in Fig. ^| the magnetocrys- 
talline anisotropy energy £w s (A so ) for the bandfillings 
nd=l and n<2=2.5. In Fig. |3|, however, we see a discon- 
tinuous change of the slope for nd—2.5, because a level 
crossing (see Fig. ^) has occurred between the levels 5 
and 6 for A so =0.04 eV. The inversion of these levels is 
indicated by the labeling in Fig. ^, for A=0.05 eV. Thus 
E an is(nd) can change its sign because of the level cross- 
ings. To explain for 71^=1 the dependence of E anis on 
SOC, we have to study the contribution of level shifts 
for different magnetization directions. We find that the 
shifting of nondegenerate states is for both magneti- 
zation directions proportional to \ 2 so . For the case of 
an perpendicular easy axis, the prefactor is larger. The 
shifting of degenerate states for magnetization parallel 
to the bond axis is in all cases proportional to A so and 
for magnetization perpendicular to the bond axis pro- 
portional to \ 2 so . To illustrate this we show in the inset 
of Fig. H the dependence of the level shifts on A so for 
the nondegenerate level 1 (this corresponds to n^=0.5) 
and the degenerate level 2 (n.d=l), for both magnetiza- 
tion directions (8 = 0;7r/2) for the diatomic pair. For 
8 = 7r/2 and small SOC there is a flat slope of the energy 
of both levels, while the slope increases with increasing 
X so . Otherwise, the slope of the energy of both levels for 
8 = is increasing less strongly with increasing X so . Thus 
in Fig. H the magnetic anisotropy energy E anis (nd = 1) 
decreases for larger SOC with increasing A so , as a re- 
sult of the different position of the energy levels for the 
two magnetization directions entering the calculation of 
E a nis{nd)- This yields another reason for the sign change 
of E anis (n d ). 

An other important consequence of the levelcrossing 
can be seen in Fig. [I]. The magnetic anisotropy energy 
E a nis{nd) of the quarter-band filling (rid=2.b and 7.5) 
of both subbands depends sensitively on the hopping 
parameters, because the level crossings then occur for 
smaller values of A so (decreasing t) or, in the other case 
of increasing t for larger values of X so . 

In Fig. [| (dashed curve), E an i S (rid) is shown as a func- 
tion of 3<i-band filling nd for the five-atom Fe ring. In 
the five-atom ring we find a translational symmetry of 
Eanis{ n d) with respect to half filling. This symmetry 
behavior of E an i S (rid) is found in all rings with a small 
odd number iV=(3,7,9,ll,..) of adatoms. In rings, with 
an even number of atoms and all chains irrespective on 



whether N is odd or even, we find a mirror symmetry 
of E an i S (nd) with respect to half filling. This behavior is 
also valid in the tight-binding model and will be discussed 
again in Sec. III.B. If the number of atoms continues in- 
creasing (iV=45,47,..), we then find a mirror symmetry 
of E an i s (rid) with respect to half filling also in odd rings. 



B. Results of the tight-binding model 

The tight-binding model is better than the diatomic- 
pair model, because it additionally allows to take into 
account the following points : (i) The change of the bind- 
ing angles between nearest neighbors of adatoms in rings 
is considered in a realistic way. (ii) Due to the extra in- 
dependence of E an i S (nd) in the tight-binding model, it is 
possible to choose any magnetic-moment direction in the 
plane of the ring or chain, (iii) Furthermore, in the tight- 
binding model it is possible to couple the free-standing 
chains and rings to the substrate. This extension of the 
theory, however, has not been performed in this work. 

For all chains it is easily checked that E an i S (nd) is equal 
in both models. There is a general trend that E an i S {N) 
decreases with increasing number N of Fe atoms in chains 
(compare Figs. 1 and 8). For rings this is not the case 
which will be discussed later. 

The solid curve in Fig. |^ shows E an i S (nd) for the tight- 
binding five ring. Although this curve is different from 
the result in the diatomic pair-model, the translational 
symmetry survives. A detailed analysis yields that the 
tight-binding model has the same symmetry behavior for 
small odd rings as the diatomic-pair model. The trans- 
lational symmetry can be traced back to the absence 
of a mirror symmetry between the levels in each of the 
two spin subbands with respect to ±J ex /2: In odd rings 
and chains both spin-subbands contain an odd number 
of electron states. If there the degeneracies of both sub- 
bands resulting from the exchange splitting are all lifted 
due to the hybridization, as is the case for odd rings, it 
is not possible to find a mirror symmetry in each of the 
two subbands. Since both spin-subbands have a transla- 
tional symmetry with respect to the zero of the energy 
scale, E an i S (nd) has in this case also a translational sym- 
metry with respect to half bandfilling. In all odd chains 
always five degeneracies (angular momentum algebra of 
3(i-states) are not lifted by hybridization. The remaining 
even number of the lifted degeneracies has a mirror sym- 
metry around ±J ex /2 in both subbands, thus E an i S (rid) 
has also a mirror symmetry with respect to half band- 
filling. In larger odd rings the translational symmetry of 
E a nis{ n d) is lost. Due to the increased number of levels 
in both subbands the additional lifting of the five degen- 
eracies, by closing of the corresponding odd chain to the 
ring, will not immediately destroy the mirror symmetry 
in both spin-subbands, as is the case for small odd QC's. 

A very important point for the theory of small particles 
is the existence of shell structure. Here we discuss in par- 
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ticular the oscillations of E an i S as a function of bandfilling 
rid for rings with an even or odd number of adatoms. If 
the magnetic moment is along the layer normal, it lies 
in planes perpendicular to all bonds of the ring. But 
if the magnetic moment lies in the ring plane (assum- 
ing exchange-enforced parallel moment alignment), it is 
twice parallel to a bond axis in rings with an even num- 
ber of N and once in rings with an odd number of N. 
This is schematically illustrated in Fig. |B| for the trimer 
(5(a)), tetramer (5(b)), and hexamer ring 5(c)). In even 
rings, where TV is a multiple of four the moment lies also 
twice perpendicular to a bond axis (see Fig. 5(b)). 

We investigate exemplarily for the four-Fe atom ring 
three different configurations of the magnetic moment di- 
rection, which are shown in Fig. g[ Since in larger rings 
the exchange interaction could possibly not enforce par- 
allel alignment of all magnetic moments as assumed pre- 
viouslytj, we calculated E an i S (n d ) for parallel magnetic 
moment alignment (p), radial moment configurations 
along the bonds (i) or half of the nearest-neighbor bond 
angle (h). Since there is a (^-dependence of E an i S (nd) in 
the tight-binding model, we find the following in-plane 
symmetry in rings for the first two configurations (p) 
and (i): 

E anis (9;<P = 0;n d )=E anis (9; ( l) = 3G0/N;n d ). (14) 

Note this symmetry is absent for the case (h). In Fig. ||, 
E a nis{nd) is shown as a function of 3<i-bandfilling rid of 
the four atom Fe ring, for the three configurations (p), 
(i), and (h). Especially for the 7r-bond, the magnetic 
anisotropy energy E an i S (rid) has the lowest values for 
the configuration (p). However, this is not the ground 
state configuration for Fe with atomic bandfilling rid~6. 
For all three configurations we get a parallel easy axis 
of E an i S (Fe), while its value is the largest where the 
magnetic moments are arranged in the configuration (h). 
This is the ground state of the three configurations for 
Fe, with the lowest total energy. In order to determine 
the ground-state configuration for the case of bandfilling 
around nd=7.5, we have to notice that we get a perpen- 
dicular easy axis for the two radial configurations and 
an in-plane easy axis for the configuration (p). In this 
case both radial configurations have the same total en- 
ergy, which has to be compared with the total energy 
that results from the magnetization in-plane (p) . A close 
inspection for all discrete bandfillings between rid=l and 
rid— 8 shows that the configuration (p) has the lowest to- 
tal energy in this range. 

Within our theoretical framework we are able to ex- 
plain a recently discovered phenomenon, viz. oscilla- 
tions of E an i S (n^i during film growth with a period of 
one monolayer L3. The in-plane lattice spacing during 
epitaxial growth of Co on a Cu(OOl) single- crystal sub- 
strate is found to oscillate as a function of coverage and 
produces anisotropy-oscillationsEfl. These oscillations are 
schematically illustrated in Fig. |^, which shows the con- 
nection between the change of the lattice spacing from R 



to i?i due to the growth and the corresponding change 
of E an i s . After starting the deposition, Co islands nu- 
cleate on the Cu surface. The Co atoms are relaxing in 
the direction of the center of the island. This gives rise 
to a reduced in-plane lattice spacing. If the growth of 
the N+l-layer is completed, £' aTl i S (N-|-l-layer) nearly re- 
turns to the value of E an i S (N- layer) . For uncompleted 
layers, however, the in-plane lattice constant is very dif- 
ferent from that of complete layers and therefore oscil- 
lations occur in E anis (N-layer) as expected. As a model 
assumption of this situation the four atom Fe-ring can be 
considered. We scale all d-electron hopping parameters 
belonging to one bond of the four-atom Fe ring with a pa- 
rameter t' ', while leaving the others unchanged. Then we 
find also anisotropy-oscillations, due to the changed posi- 
tion of some levels in the energy level scheme of the four 
atom Fe ring. For decreasing t' for one bond, E an i S (n d ) 
is increasing or decreasing dependent on nd- 

Finally we would like to discuss in this paper some 
results for larger QC's and QWs, in order to reach the 
experimental situation. In Fig. |8|, E an i S (rid) is shown as a 
function of 3d-bandfilling nd of the 48 atom Fe ring (solid 
curve), for the p-configuration of Fig. |(| A comparison 
of the three configurations (p, i, h) shows only marginal 
deviations of the total energy for the 48 QC, thus we may 
assume the p-configuration to be the ground state. For 
an effective bandfilling of n e d =6.6 we find an in-plane 
easy axis with E an i S (Fe)—lA5 meV, assuming implicitly 
a hybridization of the (^-electrons with s-electrons for this 
larger system. The shape of the curve is very similar to 
that of the tight-binding calculations for free-standing 
Fe-monolayeral3E3. As can be seen from the comparison 
of Figs. ^, ||, an d |, E anis (N) is essentially independent 
of N. 

Since metal atoms aggregating during film growth 
preferably on step edges (see ref.[20] for Fe/W(110)) may 
cause one-dimensional metallization and new magnetic 
properties, it is of interest to calculate the magnetocrys- 
tallinc anisotropy of larger QWs. The origin of mag- 
netic step anisotropics results from the modified local 
symmetry of step atoms compared to surface atoms. In 
Fig. |8|, E an i s (nd) is shown as a function of 3cZ-bandfilling 
nd for the 48-atom Fe chain, for two different scalings of 
the hopping parameters (the dashed curve refers to the 
unsealed case t=l). The value of E an i S (n d ) is consider- 
ably lower compared to the diatomic pair. However, it 
is still much larger then the values for freestanding Fc- 
monolayers. The shape of the curve for t=l is very simi- 
lar to the diatomic pair (see Fig. |l|), irrespective of atomic 
number N . Thus we expect very large contributions to 
Eanis(n d ) in all chains from SOC-induced lifting of de- 
generacies. Comparison of Figs. [I] and || indicates that, 
upon decreasing t, the shape of the curve will change 
similarly for the 48 atom QW and for the dimer due to 
levelcrossings. 

In order to show how sensitive the the magnetocrys- 
talline anisotropy energy depends on spin configura- 
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tion, we calculate E an i S (nd) for two 48-atom Fe chains 
which are antiferromagnetically (AF) or ferromagneti- 
cally (FM) coupled to each other (Fig. ^). In these 96 
atom systems E an i S (N) has further decreased in both 
cases. In AF coupled chains E an i S (rid) is less dependent 
on the scaling of the interchain hopping than in FM cou- 
pled chains. Note that the shape of E an is(nd) for AF 
coupled chains is similar to that of the 48 atom Fe chain. 
In this case the shape will not change if the interchain 
hopping tl is varied, in contrast to FM coupled chains. 
Thus the breakdown of the one dimensionality in two 
FM coupled chains will change E an i S (rid} very strongly. 
We find in this case a perpendicular easy axis with 
E a nis(Fe)=~ 0.36 meV (tl=l). This value is of the same 
order of majpwtudc as for Fe-monolayers. Experimental 
results showEil that surface and step anisotropies are of 
the same order of magnitude, both exhibiting symmetry- 
breakingj as was introduced by Neel's phenomenological 
modeE3'El This is in agreement with our results for two 
FM coupled chains. The opposite sign and same order 
of magnitude for surface and step anisotropies is repro- 
duced. 



IV. SUMMARY AND OUTLOOK 

A calculation of the magnetocrystalline anisotropy 
energy E an i S (rid) for free-standing chains and rings 
of Fe adatoms is performed. SOC is included non- 
perturbatively. E an i s (N) for ID chains is larger then 
for 2D rings, but in both cases E an i S is enhanced co- 
pared to Fe-monolayers. Generally, very large contribu- 
tions to E an i s (rid) result from the SOC-induced lifting 
of degeneracies. The shape of E anis (N) is in all chains 
very similar. In rings E an i S (N) is essentially indepen- 
dent of N, while in chains E an i S (N) generally decreases 
upon increasing N . For all chains and even rings, there 
is a mirror symmetry of E an i s (rid), while for small odd 
rings we find a translational symmetry. In small rings, 
there are odd-even oscillations of E an i S (N), which might 
be of technological relevance for magneto-optical switch- 
ing. Larger rings approach the results for freestanding 
Fe-monolayers. E an i S (nd) is found to be more sensi- 
tive to interchain hopping variation for two ferromagnet- 
ically than for two antiferromagnetically coupled chains. 
A generalization of Hund's third rule on SOC is found. 
Oscillations of the magnetocrystalline anisotropy energy 
due the growth of Co on Cu have been explained. 

Our results indicate that low-temperature scanning 
tunneling microscopy is a promising new tool for the 
preparation of atomic-scale nanostructures. These 
nanostructures exhibit interesting magnetocrystalline 
anisotropy properties which can be tailored by choos- 
ing appropriate atomic species and atom number, ring 
or chain conformations, and a suitable substrate. It will 
be of theoretical interest to incorporate the substrate de- 
grees of freedom in the investigation of the magnetocrys- 



talline anisotropy of quantum wires and corrals. 

It is then to be expected to establish many of the 
phenomena successfully demonstrated in thin film mag- 
netism, also for one-dimensional magnetic systems, such 
as enhanced moments, giant magnetoresistance, quan- 
tum well oscillations, giant nonlinear Kerr rotation, or 
magnetic reorientation. 
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APPENDIX: HUND'S THIRD RULE. 

The curve of E™ i ipar (rid) changes it sign at rid=2.5 
and 7.5 (see Fig.l). This sign change occurs exactly, 
when in each rase the two subbands are half filled. 

Wang et alB also found this behavior for SOC between 
the opposite spin states for the dimer which can be under- 
stood from the following equation (z denoting the sym- 
metry axis and x an arbitrary direction in the perpen- 
dicular plane), 

AE ud = E ud (x) - E ud (z) 

o~ 

ud (o~) represents SOC between the opposite spin states 
(represents occupied spin-down states), and AE ex is the 
exchange splitting. The sign of AE u d depends only on 
the axial component of the angular momentum. Three 
states (L z = 0, ±1), which get first occupied, give negativ 
contributions and two (L z — ±2) contribute positively. 
Due to the hybridization, the two antibonding states with 
L z = ±2, which contribute positively, will be occupied 
next, and therfore it comes to a sign change. 

This-|behavior is also valid in larger chains and mono- 
layersli3E2l. If one spin-subband is half filled, there is L=0 
and from this results no contribution to the anisotropy 
energy. The antibonding states will be occupied in re- 
verse order and thus there is a sign change at nd=2.5 
and 7.5. 
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FIG. 1. Dependence of the magneto-crystalline anisotropy 
energy E anis on the 3d-bandfilling n d for a diatomic pair. 
While the dashed curve refers to the hoppings given in the 
text the dotted curve refers to hoppings scaled with an over- 
all factor of t=0.7. Negative values of E anis (n d ) correspond 
to a perpendicular easy axis. The solid curve show the contri- 
butions of E a l ™ll par (n d ) to E anis (n d ) from antiparallel spins. 
As can be seen the main contribution to E an is(n d ) results 
from E p ^ r ig (n d ), due to the large exchange splitting J ex in Fe. 
The zeros in E a l ™ t i l par (n d ) at n d =2.5 and 7.5 indicate a sign 
change inposed by a generalization of Hund's third rule (see 
text). 



FIG. 8. Dependence of the magnetocrystalline anisotropy 
energy E ania on the 3d-bandfilling n d for a 48 atom ring (solid 
curve) for the tight-binding model, with the configuration 
that all moments in plane are parallel to each other. The 
dashed (dotted) curves shows the dependence of E an i s (n d ) 
for a 48 atom chain for two different scalings of the hopping 
t=l (t=l/3). 

FIG. 9. Dependence of the magnetocrystalline anisotropy 
energy E anis on the 3d-bandfilling n d for two antiferromagnet- 
ically (dashed curve) and ferromagnetically (solid and dotted 
curve) coupled 48-atom Fe chains for the tight-binding model. 



FIG. 2. Energy levels of the diatomic pair for the two mag- 
netization directions E m (9 = 0) (magnetization in bond di- 
rection) and E m (9 — n/2) (perpendicular to the bond axis). 
For simplicity only the spin-up states of the energy scale are 
shown in the plot. For magnetization perpendicular to the 
bond axis the degeneracies are very weakly lifted by SOC, in 
contrast to the case of magnetization in the bond direction, 
where SOC lifts the degeneracies so strongly, that level 5 and 
6 intersect. 

FIG. 3. The dependence of the magnetocrystalline 
anisotropy energy E anis on \so (eV) for n d =2.5 (solid line) 
and 1 (dashed line) for the diatomic pair. For n d =2.5 a dis- 
continuous change of the slope occurs due to the crossing of 
levels 5 and 6 in the case of a magnetization parallel to the 
bond axis. For n d =l E anis decreases with increasing \ so . 
This behavior is explained in the inset of Fig. 3 (see text). 
The inset shows the shift of levels 1 (n d =0.5) and 2 (n d =l) as 
a function of \so for the diatomic pair and magnetization in 
bond direction (9 = 0) and perpendicular to bond (6 = tt/2). 



TABLE I. The nonvanishing energy integrals in terms of 
two-center integrals, which we used for the tight-binding 
model. 



Ey 

E z 
E x 
E x 
E, 
E,. 



x y , x y 4 



Zl 2 m 2 {V dda ) + (I 2 + m 2 - 4l 2 m 2 ){V dd «) + (n 2 + l 2 m 
3m 2 n 2 (V dda ) + (n 2 + m 2 - 4m 2 n 2 )(V dd7I ) + (I 2 + m 2 
3n 2 l 2 (V ddrT ) + {n 2 +l 2 - 4l 2 n 2 ){V dd7r ) + (m 2 + l 2 n 2 
\lm{l 2 - m 2 ){V ddtT ) + 2lm(m 2 - l 2 )(V dd ^) + \lm{\ 2 - 
V3lm[n 2 - \{l 2 + m 2 )](V dda ) - 2V3lmn 2 (V dd7T ) + ±VSml 
(I 2 - m 2 ) 2 (V dda ) + [I 2 +m 2 - {I 2 - m 2 ) 2 ](y dd ^) + [n 2 + ±( 
±V3(l 2 ~ m 2 )[n 2 - \{l 2 + m 2 )](V dda ) + i^3(l + n 2 )(l 2 
[n 2 - \(l 2 + m 2 )] 2 (V dda ) + 3n 2 (l 2 + m 2 )(V dd7I ) + |(/ 2 + 



FIG. 4. The dependence of the magnetocrystalline 
anisotropy energy E an i s on the 3d-bandfilling n d for the five 
atom ring calculated in the diatomic pair model (dashed 
curve) and the tight-binding model (solid curve). In both 
models, there is a translational symmetry with respect to 
half filling n d — 5. In the tight-binding model all moments 
in-plane are choosen to be parallel to each other for this case. 



FIG. 5. Odd-even oscillations of the rings (schematic). 



FIG. 6. Dependence of the magnetocrystalline anisotropy 
energy E an i 3 on the 3d-bandfilling n d for a four atom ring 
calculated in the tight-binding model for the three displayed 
magnetic moment configurations. In (p) all moments are par- 
allel to each other, while in (i) and (h) the two radial config- 
urations are given, where all moments are in the bond direc- 
tions or half to these. 



FIG. 7. Schematic sketch of anisotropy oscillations during 
the growth of Co on Cu. We expect that the change of the lat- 
tice spacing from 7? to Ri due the growth of the (N+l)-layer 
will produce oscillations of E anis , which are found in the cal- 
culation for the four atom Fe ring, used to simulate this sys- 
tem. 
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